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Abstract 

In this paper we study the genealogical structure of a Galton-Watson process with neutral 
mutations, where the initial population is large and mutation rate is small [3]. Namely, we 
extend in two directions the results obtained in Bertoin’s work. In the critical case, we 
construct the version of Bertoin’s model conditioned not to be extinct, and in the case with 
finite variance we show convergence after normalization, of allelic sub-populations towards a 
tree indexed CSBP with immigration. Besides, we establish the version of the limit theorems 
in [3], been for the unconditioned process and for the process conditioned to non-extinction, 
in the case where the reproduction law has infinite variance and it is in the domain of 
attraction of an a-stable distribution. 

Keywords. Branching process; Nentral mntations; Q-processes; Regnlar variation; Domain of 
attraction of a-stable laws. 


1 Introduction models and main results 

A Galton-Watson process models a population where at every generation each individual repro¬ 
duces according to the same distribution, independently of the others and then dies. A number 
of variants, involving different types of conditioning and limit theorems, are core of branching 
processes theory. For instance, when the process dies with probability one, Yaglom (1947) proved 
that the distribution of the process conditioned to non-extinction exists, under some assumptions 
on the moments of the reproduction law. The proof was simplihed and the moment assumptions 
removed by Joffe (1967), Seneta and Vere-Jones (1967). More generally, Lamperti and Ney (1968) 
introduced the Q-process, conditioning on not being extinct in the distant future and being ex¬ 
tinct in the even more distant future. For background about branching processes we refer to [1], 
[12] and US). 

As an extension of the Galton-Watson model, Bertoin [3] studied the so called Galton-Watson 
process with neutral mutations; this emerges assuming that the mutations modify the genotype 
of individuals but not the dynamic of a population modeled by a standard Galton-Watson. Since 
mutations appear in the ancestral lines of the population, each individual begets children that 
do not necessarily inherit its genetic type (allele). In addition, we suppose that the population 
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has infinity alleles, that is, each mntation event originates a different allele. We denote the size 
of a typical family size by := where are non-negative random variables 

which determine the number of clones and mutants children of a typical individual, respectively. 
We exclude the degenerate cases = 0 or = 0. 

In [3], Bertoin established asymptotic features on the genealogy of allelic sub-families in a 
Galton-Watson process with neutral mutations. In his development, the genealogy of the popu¬ 
lation is described by a planar rooted tree where the mutations are represented by marks in the 
edges between parents and mutant childrens. The vertices with n marks in their ancestral line 
are associated with the called n-type individuals; those individuals with parent of (n — l)-type 
are known mutants of the n-type. We denote by the total population of individuals of the n-th 
type and by the total number of mutants of n-th type, with the convention that mutants of 
the 0-th type are the ancestors, that is Mq = a, P^-c.s. Observe that for every n each mutant 
of the n-type begets a sub-family, which is independent of the others and has the same behavior 
that the original tree. In fact, this is a consequence of the general branching property, which 
states that conditionally on the set of children of a stopping line, the families that those beget 
are independent copies of the initial tree. The concept of stopping line was introduced formally 
by Chauvin [S]. Namely, a line is a family of edges such that every branch from the root contains 
at most one edge in that family. One may think that a stopping line is a random line such that 
the event, “an edge is in the line”, only depends on the marks found on their ancestral line. An 
important consequence from the general branching property is the following lemma. 

Lemma 1.1 (Bertoin [3], Lemma 1). Under 


{Mn : n e Z+}, 


is a Galton-Watson process with reproduction law Pi(Mi G ■). More generally, 


{{Tn, Mn+i) : n G Z+}, 


is a Markov chain, with transition probabilities 

Va{Tn = k,Mn+i = I \ T„_i = i, = j) = Pj(T q = k,Mi = 1), i,j,k,l G Z+ and j < k. (1.1) 

Remark 1.2. Since the mutants of the n-th type are also individuals of the n-th type, the 
transition probability in (1.1) is zero when j > k. 

For notational convenience (^i) ■ hjjk,l G Z+} denotes the n-step transition probabili¬ 

ties of {(T„, M„+i), n G Z+}, that is 




Pa(Zr-)-n k, l\Tr hj), i,j,k,l G Z_|_, n G N. 


( 1 . 2 ) 


Remark 1.3. For i,j, k, I G Z+, P^ij) (j. depends only on the second coordinate, actually is not 
difficult to prove using induction that 


PUiKO = E PSL.)P:^-.(ro = = ;), 

jn — 1 — 1 


(1,3) 


where jo = j and {Ppj) ■ hJ ^ ^+} denotes the n-step transition probabilities of {M„,n G Z+j. 
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Let U be the finite sequence of integers 

U := IJ 

A:eZ+ 

where N = {1,2,...} and = {0}. We consider a locally finite rooted tree A = {Au : u G U} 
called tree of alleles. According to Ulam-Harris-Neveu labeling the root is 0, a vertex at its level 
A; > 0 is M = (til, ...,Uk) and uj = {ui, ...,Uk,j) is the j-th children. We denote by |m| the level 
of the vertex u. A tree A, is constructed recursively in |3] taking A 0 = Tq and Auj as the size 
of the j-th allelic sub-population of the type |m| -|- 1 which descend from the allelic sub-family 
indexed by the vertex u, with the convention that in the case of ties, sub-families are ordered 
uniformly at random. Since the transition probabilities of the chain {(T„, Mn+i) : n G Z} depend 
only on the second coordinate, the tree of alleles enjoy a kind of branching property, the formal 
statement is given in the following result, where the notation [du },) means that the d^-tuple has 
been rearranged in the decreasing order of the first coordinate, by convention, in the case of ties, 
the coordinates are ranked uniformly at random. 

Lemma 1.4 (Bertoin [3], Lemma 2). For every integer a > 1 and k > 0, under eonditionally 
on {{Au, du) : \u\ < k), for each vertex u at level k with Au > 0, the family of variables {{Auj, duj) '■ 
1 < J < du) o,re independent with distribution (Tq, under Pi. 

It is important to observe that if du '■= max{j > 1 : Auj > 0}, agreeing that max0 = 0, then 

Tk= Au and ^ du- (1.4) 

\u\=k \u\=k 

Hence given a population with neutral mutations, the tree of alleles records the genealogy of 
allelic sub-families together with their sizes. Also, the size of their generations is a Galton-Watson 
process. 

Besides, we consider for every integer n, a Galton-Watson process started from n ancestors, 
with a fixed reproduction law, which is critical with finite variance and mutations occur at random 
with rate 1/n. The main result of [3] states that as n tends to infinity, the rescaled tree of alleles 
converges in the sense of finite dimensional distributions towards a process {Zu : u G U}, 
called tree-indexed GSBP with values in (0, 00 ), such that the vertices on a level are distributed 
as the family of the atoms of a Poisson random measure conditionally on the preceding level. 
This result gives another point of view to the convergence of a Galton-Watson process rescaled 
towards a discrete time continuous state branching process, GSBP for short. 

A first goal in this paper is to construct the chain {(T„,M„+i) : n G Z+} conditioned on 
non-extinction of the mutants, so we are interested in the situation where 


T = inf{?7, > 1 : = 0} < 00 , Pq-c.s. 


(1.5) 


That is the purpose of the following Theorem. 


Theorem 1.5. Let a G Z+ and Tn the natural filtration of the process {(T„_i,M„) : n G N}. 
There exits a probability measure P^ that can be expressed as a h-transform of Pa using the {J^n)- 
martingale 


Yn = 


Mnq 


Mn — a 




4n<T}, 
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where f{y) = Ei{y^^) and q = P(0 < T < oo). 

That is 

= neEl. 

a 

Furthermore, is the law of a Markov chain {(T^, n G Z+} with n-step transition prob¬ 

abilities, 

Q'(i,j),{k,i) = jj\q))n^irj),(k,ip (1-6) 

where {PJl i) '■ i,j,k,l G Z+} denotes the n-step transition probabilities 0 /{(T„, M„+i), n G 
Z+}. 

Next we ensure that the process dehned in the above Theorem is distributed as {(T„, Mn+i) : 
n G Z} conditionally on non-extinction of mutants. Actually, thanks to a relation that we establish 
in Section]^ between the generating functions of (Tjj, M„+i) and M„, we can use classical methods 
and prove the following result. 

Theorem 1.6. Suppose that < 1 and < 1. 

i) The Yaglom limit 

lim P(T„_i = i, Mn = j\n < T < 00 ), 

n^oo 

exists and has a generating function (p[x, y) such that for all n G N, 

m^tp{x,y) = f{(pn{x,y)) - f{(pnix,0)), x,|/G[0, 1]. (1.7) 

a) Let a G Z+ and n fixed. The conditional laws of the process {{Tk, Mk+i) : 0 < k < n — 1} 
under Pa(- \n + k <T < 00 ) converge towards the probability measure P^, in the sense that 
for any n > 0 

lim Pa(A| n + k <T < 00 ) = Pa(A), VA G JTi- (1-8) 

fc^OO 


Remark 1.7. In the situation < 1, according to Corollary 1 of [5], the hypothesis E(,^A)) < 

1 is equivalent to the condition (1.5). 


Furthermore, in Section |3.3| we will see that just as in the case of Galton-Watson processes 
conditioned to non-extinction, the dynamic of {(Z([, : n G Z} admits a description in terms 

of immigration of mutants of the n-th type into the population. 

Besides we investigate, in the same sense of [3] but on a complimentary class of reproduction 
laws, the asymptotic behavior of the population. Namely, we consider a Galton-Watson process 
{Zif,n G Z+} such that the reproduction law 


7r+ = P(^(+) = 


+ ) 


is critical but, instead of assumming that it has hnite variance as in Bertoin’s paper, we suppose 
that there exists a G (1,2) such that, 

^+(i) := P(e^+^ > J) e (1.9) 
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where -RV^“ denotes the class of functions which are regularly varying with index —a at oo, see 
Chapter I in |1] for background. Note that the case where a G (0,1) is not allowed here because 
it contradicts the assumption 7r+ is critical. 

In order to extended the main result of [3] to our setting, we use that there exists a regularly 
varying function r with index 1/a such that 


r(n)P(^+ > ny) -^ c^iy 'iy > 0, 

n^oo 


( 1 . 10 ) 


where Cq, = l/F(3 — a). A proof of this fact is given in Proposition 4.1 Also, we assume that 
each child is a clone of her mother with fixed probability 1 — p and a mutant with probability p, 
so the joint law of denoted by tt = : k,l G Z+), 


= = = ^,^> 0 , ( 1 . 11 ) 

satishes 

vTfc,; = -p)V, kJeZ+. (1.12) 

In the paper [3], it has been assumed that the number of ancestors and mutation rate have the 
following behavior 

a{n) ~ nx and pin) ~ cn“^, as n —>■ oo. (1-13) 

where c, x are some positive constants. Here we will instead assume that following behaviour 


a{n) ~ xr{n)p{n) and p{n) ~ cn as n —)■ oo. 


(1.14) 


In the remainder of this paper the relation f ~ g refers to lima;_,.oo f{x)/g{x) = 1. The notation 
refers to convergence in distribution as n —)■ oo and £ j to the distribution of the 

process under The result below extends to our setting the main result in [3]. 

Theorem 1.8. If (1.9) and (l.lf) holds. Then, the following convergence holds in the sense of 
finite dimensional distributions 

jC ((((r(n))-M!.">, (r(n)p(n))-‘£i(")) : t. e U), ^ : « e U), 

where {Z^ : u E V} is a tree-indexed CSBP with reproduction measure 

u^{dy) = c'^x~^~^^^dy, ?/> 0, a G (1,2), (1.15) 


where c'^ = a ^/r(l — a ^). 

Finally we establish the convergence of the hnite dimensional distributions of the rescaled 
chain {T„, Mn+i}, conditioned to non-extinction of mutants, towards a continuous state branching 
process with immigration in discrete time. 

Theorem 1.9. If the reproduction law is critical, there exist sequences bi{n) and b 2 {n) such that 
the following joint convergence in the sense of finite dimensional distributions hold: 

C ({hi{n)Tu-iM{n)Mu) : k G Z+),P^jOf) ^ ((l/c/S^) : A; G Z+), 

where {Y^ : /c G Z+j is a CSBP with immigration such that 
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a) when the reproduction law has finite variance and (1.13) holds, then its reproduction 
measure is 

c 

: exp 


u{dy) = 


\j2na‘^y‘^ 


c^y 


2(t2 


dy, 2 / > 0, 


(1.16) 


then the immigration measure is zv^dz) and (3 = c; moreover hfin) = n ^ and h 2 {n) = n 


b) when the assumptions (1.9) and (1.14) hold, the reproduction measure is ^^{dz) as defined 


in (1.15), the immigration measure is ziy°‘{dz) and fi = 1; the normalizing constants are 
given by bfin) = (r(n)p(n))“^ and b 2 {n) = (r(n))' 




The remainder of the paper is organized as follows. In Section we provide some preliminary 
facts. Section]^ is devoted to construct and to interpret the process {(T„,M„+i) : n G Z} con¬ 
ditioned to non extinction. The last section is divided in three parts; the first one correspond to 
the framework where the reproduction law is in the domain of attraction of an a-stable distribu¬ 


tion. Section 4.2 includes the proof of our main result in this setting. Finally we prove Theorem 


1.9, which is the result that explains the asymptotic behaviour of the process conditioned to 
non-extinction. 


2 Preliminaries 

In this section we obtain some useful formulas for the generating function of {T^, M„+i), denoted 
for n G Z by 

(fn{x,y) := x,y ^ [0,1], 

where for notational convenience (pi{x,y) := (p{x,y). Observe that the generating function of 

is 

fniy) ■= g^n{l,y), 2 /g[ 0 , l], (2.1) 

and as before we denote fi{y) =■ f{y). 

According to classical theory of branching processes, the extinction probability of the Galton- 
Watson process {M„,n G Z+}, that we denote by q, is the smallest root of f{y) = y, which is less 
or equal than one depending on whether, the mean of the reproduction law, m := E(Mi) is > 1, 
or < 1, respectively. In order to avoid trivial cases, we assume throughout that 

HI) P(Mi = 1) > 0, 

H2) P(Mi = 0) -|- P(Mi = 1) < 1, and P(Mi = j) ^ 1, for any j. 

We also know that the n-step transition probabilities {Pp j) '■ hj ^ 2+} of G Z+j satisfy 

OO 

i>l. ( 2 . 2 ) 

i=o 

Now for a Galton-Watson process with neutral mutations, let g be the generating function of the 
reproduction law of a typical individual, that is 

g{x,y) := E(a;^^"V^'™^), x,y e [0,1]. 
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Proposition 1 of [3] ensure that the law of (Tq, Mi) can be obtained applying the Lagrange 
inversion formula to the equation 


f(x,y)=xg(f(x,y),y), x,!/€lO, 1|, 


(2.3) 


which gives the identity 


Pa(Po = k, Ml = 1) = Y^k-ah k > a>l and I > 0, 

k 


(2.4) 


where denotes the k-th convo lutio n of vr, as defined in ( 1.11 ). 

The statement in the remark 1.7 is a consequence of (2.3). We also have, < cx) if 

and only if E(M^) < oo. From ( |2.4[ ) we can write the hypothesis (HI) and (H2) in terms of the 
reproduction distribution of a typical individual. Besides, we have a similar identity to (2.2) for 
the process {(T„, M„_|_i), n G Z+} 


/c,/=0 


(2.5) 


where P” denotes the n-step transition probabilities. We get this equality by induction, for we 
apply the Chapman-Kolmogorov equation to express the (n + l)-step transition probabilities in 
terms of the transitions in one step and use (1.1). A simple but key relation for our analysis is 

( 2 . 6 ) 


'fn(x,y) = f„-iMx,y)), I, !/e[0,1). 


Due to (2.1), the proof of the identity is equivalent to establish the following equality 

y>n(x,y) ^ ip„-i(l,ip(x,y)), x,ye | 0 , 1 ], 


(2.7) 


which holds from the following standard calculations 
iPn{x,y) = \ T„_2,M„_i)) 

OO OO OO 

= ^ Pi(T„_2 = i, Mn-i = i) x^y^Ei{Tn-i = k,Mn = l\ T„_2 = i, M^-i = j) 

ij=0 k=j 1=0 

oo oo oo 

= Pi(T„_2 = M^_i =j)Y,Y 1 = k,Mi = 1) 

i,j=0 k=j 1=0 

oo 

= Pi(Tn-2 = i, Mn-i = y)y 

i,j=0 

= (pn-i{l,(p{x,y)), 

where we use that M„) : n G N} is Markovian, Lemma O and the branching property. 


3 The process conditioned to non-extinction 

This section is devoted to study the process {(T„,M„+i),n G Z+j conditioned to non-extinction. 
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3.1 Construction 


Here our aim is to prove Theorem that ensures the existence of a Markovian process that we 
shall understand as the chain {(T„, n G Z+}, conditioned to the non-extinction of mutants 

in the population. 


Proof of Theorem An application of the monotone convergence Theorem along with an ele¬ 
mentary computation shows that 

Moreover, the following identity is deduced form the branching property of the Galton-Watson 
process {M„ : n G Z+} and the properties of its generating function 

I,., = <-■/;(?). 

The latter and former identities imply in turn that 


Then by the Markov property, 

Combining the latter with the fact that /^(g) = [f'{q)]^ (see Athreya and Ney [1], Lemma 3.3), 
we have that 

^ ~ {nq)Y ’ 

is martingale. Now from the theory of h-transform there exists a Markovian process (see 0 ). 
that we denote by {(2^, n G Z+} whose law satishes 


= io, Ml = ji, • • ■ , = C_i, Ml = jn) := Pa{An) 


Jnq- 


'Jn O' 




(3.1) 


where for every n G N, and 

An {Zq f 0) M^i jly ' ' ' ) T'n—1 in—lj M^i jn\ i Yy jlj ■ • ■ y ^n—ly Jn^EJ. (3.2) 


3.2 Conditional laws 


This subsection is devoted to prove Theorem 1.6, in this aim the keystone will be the generating 
function of {(T^, Mn+i) ■ n G Z}, hence some of the results given in sectionwill be neccesary. 







Proof of Theorem 1.6. i) We will first ensure that the generating function converges. For this 
end observe that {M„ > 0} under the event {n < T < cx)}. Then for all x,y E [0,1], we 
have the identity 

Tn{x,y) := \n <T < oo) 

_ (Pn{x,y) - (Pn{x,0) 

1 - P{Mn = 0) 


Using the identity (2.6) and the fact /n(0) = P(M„ = 0), we get the expression. 

1 - /n-l(O) f fn-l{(p{x, y)) - /n-l(O) /„_i((p(x, 0)) - /n-l(O) 


Tn{x,y) = 


1 - /n(0) 


l-/n-l(0) 


l-/n-l(0) 


Now we take u = /„_i(0) and use that m = /'(I) to obtain the following limit. 

,, l-/n-l(0) ,, l-U 1 

lim - „ , , = hm 

n^oo 1 - /„(0) 

Besides, note that for each s the function 


If — f{u) m 


n 


1 - /n-i(s) 


l-/n-i(0)’ 

is a decreasing function. This implies that as n tends to infinity the expression 


fn-i{s) - /n-l(O) 

l-/n-l(0) 


= 1 


1 - fn-l{s) 

l-/n-l(0)’ 


has a limit, say 1 — /(s). We stress (see for instance Athreya and Ney pQ, Theorem 8.1) 
that /(s) is the generating function of m = {rhk : k G N}, the Yaglom distribution of 
{Mn : n G Z_|_}, which is obtained for all fc G N, as follows 

fhfc = hm P(M„ = k\n <T < oo). 


This fact implies that, the generating function / satisfies 

1 -/(/(s)) = -/(s)), sg[0,1]. 

The above calculations imply the following 

f{T{x,y)) -f{ip{x,0)) 


(3.3) 


ip{x,y) := lim ipn{x,y) = 


m 


Now we prove by induction (1.7). If n = 1, it is the just proved equality. Then suppose 
(1.7) holds for n = k. In order to get the identity for n = k + 1 note that by the induction 
hypothesis 

m''^'^lp{x,y) = m 1 - f{ipkix,0)) -m I - f{ipk{x,y)) , x,2/G[0, 1]. 

From where we deduce the claim by using first (3.3[) and then (2.6). 
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ii) Let us consider the event of the form in (3.2). Since {(T„, M^+i), n G Z+j is Markovian we 
have for An as in (3.2), that 


Va{An,n + k <T < oo) = Ea(lA„Ej^(l{Mfc>0}g^'“)), 


and similarly 


Pa(n + k <T <oo) = Ea(l{M„+fe> 0 }<? "+'“)• 


Then using (2.2), we get 


Va{An \n + k <T <Oo) = Pa(^n) 




i=i 


E 

J=1 


Theorem 7.4 in [T] establishes that the following limit holds 


lim 


fc—^oo P 


pn+fc 

= iiij‘(/'(9))‘9‘‘“‘". 


(*2j) 


Finally, thanks to hypothesis (H2) we can apply the previous result to obtain 


JnQ- 


•Jn Oj 


lim Pa(24„|n + k <T <oo) = Va{An) , , , a G N, 

a{r{q)r 


k^oo 


which hnishes the proof. 


Remark 3.1. In the previous proof we established the existence of a Yaglom limit when m < 
1, however similar arguments can be used to show that existences of a Yaglom limit in the 
supercritical case. 

3.3 Interpretation 

Motivated by the interpretation of a Galton-Watson process conditioned to non-extinction given 
in m, in the present subsection we provide a description of {(ZjJ, : n G Z} as immigration 

of mutants of the n-th type into the population. We start calculating the generating function of 
the n-step transition probabilities of the process. 

Proposition 3.2. The generating function of the n-step transition probabilities for the process 

6 Z+} 


is given by 


oo i_j n n-1 

^ qy) JJ qy)), x,ye [0,1]. (3.4) 
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Proof. Since x, y < 1 the generating function of (1.6) is infinitely differentiable, then an elemen¬ 
tary calculation and the use of (2.5) imply the following identity. 

,1-i 




yq 


k,l=l 


[f'iqW 


d 

Pi~\x,u)—(pn{x,u) 


u=qy 


To hnish the proof we apply repeatedly (|2.6) and the recursion fniv) = f{fn-i{y))- 


Taking a; = 1 in (3.4) and recalling the fact that the transition probabilities of {(T„, M„+i), n G 
Z_|_} depend only on the second coordinate, we can identify a Galton-Watson process with immi¬ 
gration uni. 

Corollary 3.3. If {M„,n G Z+j is critical or subcritical, then {M\ — l,n G Z+j is a Galton- 
Watson process with immigration [/, . 

Note that {M\,n G Z+} is the Q-process associated to the Galton-Watson process {M„,n G Z+} 
(see for instance dH or [U). The following Gorollary is analogous to Proposition 1 in [3]. 

Corollary 3.4. If {M^^n G Z+} is critical or subcritical, then the generating function of{Tl,Ml) 
is given by the eguation 

= —^g{^p{x,y),y), x,y E [0,1]. 
m ay 

Moreover, the distribution o/(Tq,M|) is given by 

Pl(d = k, Ml = ;) = PkW k>a>lyl>ll 


Proof. Taking n = 1 in the equality (3.4), 


Ei(a; 


T' M 
-^0 


h y ^ ^ 

G = —• 
m oy 


Then the hrst identity is obtained with the substitution of (2.3). To get the second one, recall 
the dehnition of given in (3.1), then use (2.4). ■ 

We can now give an interpretation to the process {(T[[, G Z+j, in terms of a tree of 

alleles with immigration = {Af : u E U), this will provide a description of the genealogical 
structure in a population conditioned to non extinction. The key elements are the tree of alleles 
A = (Au : M G U) and the above corollaries. 

We start dehning ZI 0 = Tg that is, the total number of individuals without mutations into the 
population, then according to a distribution with generating function f'/m, a random number 
of individuals of the same genetic type arrive. We enumerate the allelic sub-populations 

4* 

of the hrst type beget by Tg in decreasing order, with the convention that in the case of ties, 
sub-populations of the same size are ranked uniformly at random. Then, using Corollary |3.3| we 


choose uniformly at random one of the hrst type sub-families in the tree of alleles, removing it 
and replace it by a population of si ze T n which begets allelic subpopulation according to M|, 
where (Tg, M|) is given by Corollary 3.4, We continue with the construction by iteration, is 
the size of the j-th sub-population allelic of type |m| -l- 1 which descend from the allelic sub-family 
indexed by the vertex u, then we choose one of the sub-families of the type |m| -|- 1 to replace it 
for one of size Tg which begets allelic subpopulation according to M|. 
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4 Asymptotic behavior 


4.1 The a-stable case 

Consider a Galton-Watson process {2'+,n G Z+} with reproduction law denoted by 7r+, that 

7 r+ = P(e(^)+e(”*) = fc), fcGZ+. 


IS 


Suppose that 7r+ is critical and (1.9) holds. The mutations appear in the population according 


to (1.12) and we focus on the situation (1.14) described in the Introduction. Our goal in this 


section is to prove Theorem 1.8, for that end we start by describing the normalizing constant r(n) 
appearing therein. 


Lemma 4.1. If condition (1.9) holds, then there exists r{n) G such that 

/ \ -I-/ 7 \ d,y 

r[n)'K^[ndy) -^ , 

in the sense of vague convergence on (0, oo), where Ca = l/r(3 — a). In particular 


exp \—t (1 — e — Xy)r{n)7r^{ndy) > ->■ e 

i[0,oo) I 




The proof is an elementary application of standard results from the theory of Regular Variation 
(see e.g. |1] for background), but we include a proof in Appendixfor sake of completeness. 


4.1.1 An approximation for the reproduction law 

In order to link the asymptotic behaviour of the reproduction law of a typical individual with that 
of the joint distribution of clones and mutants, we hrst link their Laplace transform. Although 
in the present setting we use some ideas of the standard Tauberian-Abelian Theorem, we remark 
that it is not straightforward application of this theorem because in the present setting we consider 
sequences of measures indexed by the positive integers that change, unlike to the standard case 
where only the normalizing constants change. 

Lemma 4.2. For every positive integer n, let be the Laplace transform of 

under the measure Assume {A(n) ■. n > 0} is a positive seguence such that A(n) —)■ 0, as 

n —)■ cx). Then 

0n(A(n),6*) ~ (j)'^ ((1 — p{n)){l — +p(n)(l — e“®)) with n —>■ oo, > 0, (4.1) 

where 0+ is the Laplace transform of In particular, 0™, respectively the Laplace trans¬ 
form respectively satisfy 

0rW~0+(p(n)(l-e-^)), 

[0,^(A(n)) ~ (1 — p{n)){l — 6“'^^"^) )], as n ^ oo,y6 > 0. 
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Proof. According to (1.12), conditionally to = k the distribution of is Binomial with 
parameter {k,p). This fact implies the following equality in law 




^(+) 


(4.2) 


2 = 1 


where {Ui, i > 1} are i.i.d. random variables with common distribution that of an uniform random 
variable in (0,1). Therefore, 


0„(A(n), 0) = = k) [(1 - p{n))e + p{n)e ‘ 


k=0 


= (/)■'■(—log(l — (1 — (1 — p(n))e — p{n)e ^))). 


We conclude the proof using (1.13) and the elementary asymptotic estimate 

log(l - y) 


y 


y^O 


-y - 1 . 


(4.3) 


In the same way it is posible to establish the following estimate. 

Corollary 4.3. For every positive integer n, let il)n be the characteristic function of = 
under the measure pA"-), Then 

'0n(A(n),6*) ~ (/)■'■( (1 — p(n))(l — + p{n){l — e ^)), as X{n) - > 0,V6' > 0. (4.4) 

n^oo 


In particular tljff, respectively the characteristic function respectively of satisfies 


[V'(A(i!))~^+((1-?(«))(! 


as n ^ oo, A(n) —0 and for all 9 > 0 . 


Proof. Similarly to previous Lemma, using (4.2) we have 


V’n(A(n), 9) = '^ = k) [(1 - p(n))e'^^”^ + p{n)e'^] ^ 

k=0 

OO 

= ^ P(^*'’''^ = k) exp{k log(l — ((1 — p(n))(l — — p(n)(l — e'®)))} 


fc =0 


To conclude we apply the asymptotic estimate (4.3). 

We can now use the results above to give an estimation for the reproduction measure. 
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Proposition 4.4. For every positive integer n, let be the reproduction laws of 

and res pecti vely. Assume {y{n) : n > 0} is any sequence such that y{n) —)■ oo as n —)■ oo. 

In the regime (1.9), the asymptotic behaviour of the tail distribution of is given by 


7r^ \y{n)) ~ Ca7r~^ (|/(n)/E(^^'^)) , as n ^ oo. 
where (•) = cn,mn and Ca = l/r(3 — a). 


(4.5) 


(4.6) 


Proof. We prove the statement for clones, the mutants case is fully similar. First note that in the 
same way as in the proof of Lemma 4.1 (see the Appendix [A]), 

= f s7r^"‘{s)ds, X > 0, 

Jo 

is a measure on [0, oo) with Laplace transform C^cn such that 

A=£,_(A) = E (l - - A{<“>e-"S"”’) , A > 0. 

We now replace A by a sequence {A(n) : n > 0} such that A(n) —)■ 0 as n —>■ oo, 

A(n)2£^.n(A(n)) = 1 - C(A(n)) + A(n)(0^)'(A(n)). 


An estimate of the term follows from Lemma 4.2 


0n(^H) ~ 0^((1 - pW)( 1 - e ^^”^)), oo. 

In order to estimate we use the fact that for every fixed n, conditionally to = k the 
distribution of is Binomial with parameter (A;, 1 — p{n)), and we proceed as in the Lemma 
0| to get 

m'(Hn)) ~ (1 - p(n))e-"l”V+'((l - p(n))(l - n ^ «d. 

Putting both terms together we infer the estimate, 

A(n)=£,™(A(n)) ~ 1 - .#.+((1 - p(n))(l - + (1 - p(n))(l - - p(n))(l - e""'”')) 

— (1 — p{n)){l — — X{n)e~^^'^^)4>f'{{l — p{n)){l — e“'^^”^)), as n —)■ oo. 

Besides, 

(1 - p{n)){l - - A(n)e-^('^))0+'((l - p{n)){l - e'^^"))) -^ 0, 

n^oo 

and from ( |A.2 ) we have 

A(n)^£^(A(n)) = 1 - 0+(A(n)) + 0+'(A(n)), n oo, 


where is the Laplace transform of the measure p. defined in (A.l). From these last two displays 
we obtain 

A(n)^£^cn(A(n)) ~ ((1 - p{n)){l - e"^^”^))^£/,((l -p(n))(l - e"^^'"^)) + 0(A(n)^), as n -)■ oo. 
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Due to the estimate A(?7,) ~ 1 — e 


(A.3) implies 


in the limit as u —)■ cx), the approximation of given in 


Ca\{nfC l^cn (A(n)) ~ TT^ 


A(n)(l — p{n)) 


+ 0{\{nY), n —)■ cx). 


(4.7) 


Hence is remains to prove 


lim 


7r‘^"'(l/A(?7,)) 


= Cn 


n^oo (A(n))2£^cn(A(n)) 

In this aim, dehne for every y > 0 the following measure 

y^^{y/X{n)) 


(4.8) 


"^ ah ( 0 , 2 /] ■=m''x(n){y) = 


Cfj_cn (A('n)) 


Observe that 


[0,oo) 


e-^’d, ('ph(M) = ) f = ‘^r , v» > 0, 

V C^cn (A(n)) J (A(n)) (A(n)) 


from the previous display and (4.7), we get 

V0>O. 

n—>-oo 

Writing now in terms of the gamma function we get that 

1 


Q-{2-a) _ 


s(2-")-ie-®*ds. 


r (2 — a) Jo 

and since the convergence of the Laplace ransform implies the weak convergence (see 0), we have 


mTin)iy) —^ 


n->-oo r(2 — a) 




,2-a 


r(3 — a) 


(4,9) 


Besides, by the dehnition of the measure we can obtain that for any ?/ < 1, the following 
inequality 

{ll\{n)fy{l-y)n‘^-{y/\{n)) < /i-(l/A(n)) -/.-(|//A(n)) < (1/A(n))^(l - 2/)7r-(l/A(n)), 


holds. Now, due to (|4.9), we obtain 


1 — 2/^ “ . „ 7r^”(l/A(n)) 7r'^'^(l/A(n)) 

Cr.—^ -< hmmf .:: .. < hmsup v / v yy 


l-y n_^oo (A(n))2£^cn(A(n)) n^oo (A(n))2£^c„(A(n)) y{l - y) 
To conclude we make y t 1. 


1 - 

< Crv^:-T for all y < 1. 
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4.2 Proof of Theorem 1.8 


Before proving Theorem let us we describe the scheme of proof that we will follow. According 
with the construction of the alleles tree {Au ■ u G U}, given m G U such that \u\ = k > 1, for 
all j G N vertices uj, represent the size of the j-th allelic sub-populations of type k + 1 begeted 
by M, thus labels of vertices at level k + 1 determine the variable T^+i. Moreover for all fc G N, 
the total number of vertices at level k correspond to Besides a tree-indexed CSBP is related 
to the CSBP in the following way: vertices u G U at level |m| = k > 1, represent the sizes 
of the sub-populations at generation k in the CSBP which descend from the same parent at 
generation k — 1. So, a hrst step to prove the converge of alleles trees in Theorem 1.8 is to prove 
the Proposition below. 


Proposition 4.5. Assuming (1.9) and (I.I 4 ), we have 


C 


T, 


Ml 


fc +1 


fc G Z_ 


P 




((AA Mja 


) : fc G Z+), 


jinY r(n)p{n) J ^ 

where : k G Z+) is a CSBP process with reproduction measure z/" given in Theorem 


1.8 


Since the processes {(T„,M„+i) : n G Z} is Markovian with transition probabilities given in 


(1.1), to get the previous convergence we need the following Lemma. 

Lemma 4.6. For a G (1,2), let r“ he an 1/a-stable subordinator with Laplace exponent 


log E(e 


-Ar“^ 


= X 


(1 






(4.10) 


where is as defined in Theorem 1.8. Assuming (1.9), (1.13) and (I.I 4 ) 
i) the following convergences holds: 


C 


f _ Ml _^ pPO\ f a ay 

Vr(n)’r(n)p(n); ’ ain)J ^ KC, C), 


(4.11) 


a) under the measure the behaviour of the joint tail distribution of {Tq, Mi), is given by 

Tn Ml 


lim r{n)p{n)PM^ 


> s, 


r{n) ’ r{n)p{n) 
where denotes the tail function of the Levy measure z/“. 


> t ] = z/"(s A t). 


(4.12) 


Observe that the convergence in (4.11) establishes Poposition 4.5 in the case k = 0. Because 
a discrete time CSBP can be constructed from subordinators, that is, given a subordinator {r“ : 
t > 0}, the sequence {Z^ : k G Z+} dehned by iteration as follows: 


fe+i 

E 

2=1 


Zi = r: 


X-\-Z\-\ -’ 


A: G Z 


+ ) 


with the convention oi Zq = x for some £x x > 0, is a discrete time CSBP, we refer to |T2] for 
further details. 
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The tree-indexed CSBP Z can also be constructed from this subordinator, namely according 
to [3], Zoi = X and {Zj : j G N} is given by the sequence of the jumps of r“ on (0,x] ranked 
in decreasing order. The sizes of the subfamilies beget by the subfamily of size Zi , denoted by 
{Zij : j G N}, correspond to the ranked sequences of the jumps of r" on the interval {x,x + rfj 
ranked in decreasing order, {Z 2 j : j G N} to those on the interval {x + Hence, 

the convergence (4.11) will establish the convergence in Theorem 1.8 with u = 0, once the 
subordinator r“, is written in terms of their jumps, or by Levy-Ito decomposition, in terms of the 
atoms of the Poisson point process associated with it. This is the aim of Lemma below. 

Lemma 4.7. Let b{n) he a sequence of integers such that h[n) ~ br{n)p{n) for some b> 0. 

a) For every n G N, let : 1 < j < b{n)} be a sequence of independent identically distributed 

random variables with distribution ■ Defining for every n G N, := h^(„) 

and following convergence holds 

In ^ 7, (4.13) 

where j is a Poisson point process with intesity br], with r] the image of the measure z/" 
(given in Theorem 1.8) by the action of the map x h->■ {x,x). 

b) We have the following convergence, under the measure 


Tn 


Ml 


r(n) ’ r{n)p{n) 


{6(n)L 


(ai, a2,...), 


(4.14) 


where for all /c G N, a*, = (a^, Ok) with {afcjfceiN, the atoms of a Poisson random measure on 
(0, oo) with intensity bu° ranked in decreasing order; the measure z/" is given in Theorem 

Proof. From the construction : 1 < j < b{n)} is a sequence of independent random variables. 
Besides, the convergence (4.12) in Lemma 4.6 implies 7 -"'^ 0 as n —>■ 00 , uniformly in j. Then 


sup A 1|) -)■ 0. 

j 

Hence, according with the definition given in chapter 4 of [Hj, we have that { 7 ]"^ : 1 < j < b{n)} 
is a null array. Thus we will get the convergence (4.13) as an application of Theorem 16.18 of [5], 
once we verify the following condition: 

i) > 0) —)■ p{B), as n —)■ 00 , for all B E B, where p is the image of the measure 
6 z/“ by the application x h-)■ (x,x), 

ii) 'Yl,j > 1 ) —)• 0 , as n —>■ 00 , for all B G B. 

where B is the Borel cr-algebra of [0, 00 )^, B =: {B G H, 7 ( 95 ) = 0 c.s.}, with 7 the measure 
defined in the statement of this Lemma and the symbol d denotes the boundary of B. Observe 
that the class of sets B = (( 6 , 00 ) x P_|_) U (P+ x ( 6 ', cx))) is a vr-system which generates a A-system 
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and it is just B. Then by Dynkin’s Theorem, is enough to establish the conditions above for sets 
of the form B. In this setting, the condition (ii) holds because takes only the values 0 or 

1, for any j and n, the condition (ii) is true. To establish (i), observe the following equalities 

b{n) b{n) 

E >o) = t2 = 1) 

i=i j=i 

b{n) 

t=i 

= b{n)fxri{{r{n))~^s, {r{n)p{n))~^t), 

here we recall the identity pn{x^y) = > a; or Mi > y). Assuming that h{n) ~ br{n)p{n) 

for some 6 > 0, from Lemma 4T and the last equality we have 

b{n) 

V > 0) ^ bu^is A t)). 

' ^ n^oo 

3=^ 

To get the hrst convergence in the Lemma, it remains to observe that the equality 6z/"(sAt) = ri{B) 
holds. But this follows from the following equalities. 


p{dx, dy) = b 


IB 


' {x,x)&B 


u^{dx) = b 


(sAi,oo) 


z/“(da;). 


We will now prove the convergence (b). For z = 1,2, Xij^ denotes the z-th coordinate of the 
sequence that appear in the statement (a). Assuming that := (Xu\x^'') ^ •= 

if and only if x^u ^ or ^ X 2 i^ iof as dehne j\ as the index where the 
maximum integer of the sequence {xl"^^ : 1 < * < ^(a)} is reached. 

(n) 

Xp = max Xi. 

Similarly for k = 2, ...,b{n), let jk be the index of the fc-order statistic 

(n) 

xj = maxXi, 

where Jk = {1,..., b{n)} \ {ji, ...,jk-i}. Then observe that 

P(xS-?^ > Cl,xl? > C2,-.,X? > Cfc) = P(7n(C'l) > l,7n(C'2) > 2, ..., 7„(C'fc) > k), 


if Cj = (cj, Cj), Ci = (0,1) X (0,1) \ (0, Cj) X (0, Ci) and ci > • • • > c^. Taking now the limit when 
n —>■ oo in the equality below and using the convergence in (4.12), we have 

P(Xii > Ci,Xi2 > C2,...,Xjfc > Cfc) -^ P(7(C'i) > 1,7(C'2) > 2,...,7(C'fc) > k). 

n^oo 

This implies the desired convergence because 

[P(7(C'i) > l,---,7(C'fc) >k) = P(aj^ > Ci,aj2 > C2,...,aj^ > c^), 

where a^ = (a^, ak) with {afcjfcgiH, the atoms of a Poisson random measure on (0, oo) with intensity 
bu°‘ ranked in decreasing order; the measure iz" is given in Theorem 1.8 As before we used the 
indices jk to rank in decreasing order the sequence a^. ■ 
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4.2.1 Proof of Lemma 14.61 

For each n G N, let be a random vector such that their coordinates are obtained from 
and as follows 


S<"> = (a(n), 0) + ^ (d“> - 1. d""’) . e N, 

i=l 

where £ ItJ} are a pair of i.i.d. random variables with distribution 


A key result in the proof of Lemma 4.6 is the following: 


Proposition 4.8. In the regime {1-9}, (1.13) and (1.14) ihe normalized random walk defined by 

[r{n)t\ 

= {a{.n)/r{n)p{n),0) + ^ - l)/n, ^ > 0, 

i=l 

converges weakly 

(S[d„„j:«>0)^((i + X„«):t>0). 

where {Xtfi > 0} is an a-stable process with no-negative jumps with and characteristic exponent 
Ca|A|“. 


We pospone the proof of this result to the appendix and focus on the proof of Lemma 4.6 


Proof of Lemma 4-6 i) . From Lemma 3 of |3], we know that the hrst passage time below 0 for 
the centered random walk + has the same distribution as Tq. Moreover, 

on the one hand, because of the identities. 


d")(0) 


= inf {A: e Z+ : = 0} and 

i=l 

the random variables 

and (To, Mi), 

have the same distribution under On the other hand, we also have the following two 

<^O)(0) 1 


identities 


mi{k e Z+ : = 0} = inf{t > 0 : = 0} 


(n) 


r(n) r(n) 


and 


— rO)|'ni sO) 


r(^n 

and the weak convergence 




)(o)J’ r(n)p{n) 


SO)(0) 


(4.15) 


: t > O) ^ ((x + Xi,t) : t > 0) 
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Since X is a spectrally positive Levy process, Theorem 1 in Chapter VII of [2] ensures that the 
first passage time below —x for the process X 


r" = inf{t >0.Xt< —x}, x > 0. 


is a stable subordinator of parameter 1/a. We conclude from these facts that the claimed con¬ 
vergence holds as soon as 




r{n) 


■W (0) 


{r^,{x + Xt,t)\t=r^) . 


(4.16) 


But according to Theorem 13.6.5 of [H] about weak convergence of first passage times and under 
shoots and overshoots, when there is convergence in Skorohod’s topology, we get that 


Moreover since we have the joint convergence 




pS 


(n) 




: f > 0 j ((x -|- Xt, t) : t > 0) 


in the sense of Skorohod’s topology, and the second coordinate is a determinist linear function, its 
is an elementary exercise to extend the above mentioned result of to get that the convergence 


in (4.16) holds 


Proof of Lemma f.6 ii) . We will apply the same techniques used in the proof of statement (ii) 
in Lemma 4 of [3]. Let us start observing that for every x, ?/ G K 


^OO /*oo 


e^^*y = st / e ^yi{^^u,y<v}dudv, s,t> 0 . 

Jo Jo 

Thus Fubini’s Theorem implies that for any random vector (X, Y) the following identity holds. 

poo poo 

1 - 

In particular. 


= st / e *^P(X > M or V > v)dudv, > 0. 

Jo Jo 


f*oo roo 


l_[£-p(”-)(g ^ I j ^ su (r(n)p(n)) '^v)dudv, s,t>0, 


'0 JO 


where fin{x,y) := P^*'”^(To > x or Mi > y). Hence by the branching property, 

^ — st I I e~^'^~^'^fLn{{'r{n))~^u,{r{n)p{n))~^v)dudv 


oo poo \ a{n) 

---- 


'0 JO 


Now, according to (i) the left term in the previous display converges as n —)• oo towards 


(4.17) 
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Besides the right term converges as n —)■ cx) towards 


exp (-X (1 - ) , 


here we use (i) in the Lemma 4.6, together with (1.13). Therefore 


(4.18) 


Taking logarithms in (|4.17|) and (|4.18|) we obtain 

poo 

lim sta{n) / 


e *’'/i„((r(n)) (r(n)p(n)) ^v)dudv = x / (1 — e 


JO JO 

Hence it only remains to see that the line above is equal to 


xst 


^-sy-tz-a 


iy°‘{y A z)dydz. 


'0 JO 


For that end we observe the equality 


Az)dydz= / u^{du) 


^ or u>z}dydz ^ 


'0 JO 


'0 JO 


and we obtain the claimed identity by uniqueness of the Laplace transform. 


4.2.2 Proof of Proposition 

Observe it is enough to show the convergence of the Laplace transforms associated with the hnite 
dimensional distributions of each processes. So, we will briefly deduce the Laplace transform of 
a CSBP (in discrete time) {Zk : k G Z+}. From [12], we know that the transition probabilities of 
this branching process are characterized, for every k G Z+ and A, ^ > 0 by 


= z) = 

where k is the cumulant of a subordinator without drift, that is 

fi:(A) = f {1 — e~^^)id{dx), 

J (0,oo) 


where -d is the Levy measure. Applying successively the property (4.19), we obtain 

Si > 0, z = 1, 2 , ..., k. 

where by dehnition lo{s) = s and 

^n—i T 1 ('Sn—i+l)))5 ^ ^ ^+' 


(4.19) 


(4.20) 


(4.21) 
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In particular, {2’^“ : k G Z+} is a CSBP whose transition probabilities are characterized by a 
subordinator r", with cumulant given by (4.20) with dehned in Theorem 1.8 Therefore 


we will prove by induction on k the following convergence: 


-p(n) 

~a{n) 




T,_i- 


r(Ti)p(n) 




^i=l 


-)■ exp{—a;K(4-i(so + for all Sj, fj > 0, i = 1, 2,..., fc. 

^OO 

(4.22) 

where I is dehned in (4.21), taking ■d = We start observing that has the same distribution 


of r“ because of CSBP’s construction. Thus the case n = 1 follows from Lemma 4.6 (i). We next 


assume (4.22) holds for n = k and prove the convergence for n = k + 1. To this end, we use the 


Markov property of {(T„,M„+i) : n G Z} and the fact that, conditionally to = j, the pair 
(T„, has the same distribution of (To,Mi) under Pj, to obtain 


-p(n) 

~a(n) 


^/c+1 


—_ 


I G r(n)p{n) 


Mi 


>^2=1 


a{n) ' 


r(n) 


“'‘-iTk-i-i , (fe, , - , 1 , , log 

I r'frj. iTilnl r'trj. iTiln. I -p 


r(n)p(n) r(n)p(n) ° r(n)p(n) 


fc + l 


g r(n) 0 r(n)p(n) 1 j j 


Due to the assumption r{n)p{n) ~ a{n)x in hypothesis (1.14), we obtain as consequence of Lemma 


4.6 (i), that 


'k+l 


eS:! n 


y2=l 


.Ml 


■“fe-i 


- T—rTi-i - , , Ml \ n-p(^) f- t^T) - , . , , i - 

p ryn) ^ ^ r{n)p{n) *' 1 ,-s^ ' ^ ( 0 Tfrij ^ r(n)p(n; ... p r[n) 

* ^a(n) ' 


'^k — l r(n)p(n) + l 


Then using the induction hypothesis with 


s'i-i + t'i — 


Si_i +ti i < k, 
l{si-i+ti) i = k, 


(4.23) 


we get 


E 


fg r(n)T> r(n)p(n) . . . g r(n) '^^-1 (tfc+K(«fe+tfc + l )) Affc 


a(n) 


-G exp{-xK(4-i(so + f'l))}. 


This concludes the proof because of the recursive dehnition of k given in (4.21) together with the 
choice of + f', 

^(4-2(Si + ^ 2 )) ~ ^(4-i(si + ^ 2 )), 

as consequence of /fc_i(sQ + ct[) = 4(so + ^i)- 
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4.2.3 Proof of Theorem 11.81 

We will establish 


L 


Ai'n) A'n-) 

Uj 


W « C! I ^ (((^.. ^») : l«l < ^) : *= e ZJ . 


(4.24) 

Actually by the monotone class theorem, it is enough to show for non-negative measurable con¬ 
tinuous functions fi,fk 


E‘ 


■p{n) 

'a{n) 




2=1 




Y\_fi{2u,Zu : |m| < i) 


2 = 1 


where Qa, is the law of a tree-indexed CSBP started with an initial population of size x contracted 
from a subordinator {rf : t > 0}. Let us prove it by induction on k. The case fc = 1 is given in the 
convergence (4.14). Assuming the case k, we will prove the c onve rgence for A; -|- 1. Conditioning 

we have 


with respect to Tk = o'(A^\ : |m| < k) and using Lemma 


1.4 


^p(n) 

a{n) 


^k+1 


O n (r(n)p(n))-^dP : |«| = ») | 


y2=l 


= E- 


■p{n) 

a{n) 


n /*((^(^)) (r(n)p(n)) : |m| < i) 


2=1 


E?”^ (fk+i{{{r{n)) (r(u)p(n)) ’1:1^1=^) 


Besides, by induction hypothesis ~ r{n)p{n)Zu with |m| = k, therefore when u —>■ cx) in the 
previous equality we obtain 






7 

ui 


|m| < 0Ei(/fe+i((a'i,a2,...))) 


2=1 


where a). = (a)., a'^) are the atoms of a Poisson random measure on (0, oo) with intensity hZ^i'^, 
repeated according to their multiplicity and ranked in the decreasing order. Due to the dehnition 
of a CSBP tree this concludes the proof. 


4.3 The conditioned to non-extinction case 


This section is devoted to prove Theorem 1.9 Following the same strategy of Proposition 4.5 


we shall establish by induction the convergence of Laplace transforms of the hnite dimensional 
distributions associated with the processes involved. In this aim, we hrstly deduce the Laplace 
transform of the hnite dimensional distributions of a CSBP with immigration, with mechanism 
{■d, l), : n G N}, dehned for every u G N as follows 


“^n-l-l 


= ^n{Z, 
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where {/„ : n G N} is a sequence of nonnegative random variables with common probability 
measure u, which determine the distribution of individual immigrants arriving in the population. 
Let us denote its Laplace transform of uj by l, i.e. 


— I ^ ^^u{dx), A > 0. 

Jo 


(4.25) 


Let {Tn(t) : t > 0}„>o be a sequence of independent subordinators (without drift) and also 


independents of In, with the same distribution and Laplace transform given in (4.20). Thereby 


This previous equality together with the Markov property imply that 

k-l 

= for all s, > 0 ,* = 1 , 2 ,...,k. 

4=0 


Thus the proof of statement (i) of Theorem 1.9 require to establish for all Si, ti > 0, i = 1,2,..., k, 
the convergence below holds 


n 


k-l 


e 41^ 


-T.i-UVfi 


-)■ exp{-xK{lk-i{so + C^ti))} TT L{lk-i{si-i + c^u)), (4.26) 

n -t- J- 


^i=l 


i=0 


where k and / are respectively defined in (4.20) and (4.21), taking in partiular D = u given in 


(1.16). To obtain (ii) of Theorem 1.9, the previous convergence is proved with'd = z/" defined in 
(1.15). The following Lemma establishes the above convergence in the case k = 1. In their proof 


we use the reference | 3 ] to justify the statement corresponding to = u and previous results for 
the associated to d = z/". 


Lemma 4.9. If (1.13) holds, then we have the following convergence 

C (^(6i(n)To,62(n)Mi),P(([”j^j ^ ('r,c^r), 

where r is a random variable with Laplace transform where n{s), l{s) are given in 

(4-^5), according to 


a) if the reproduction law has finite variance cx^, D = c where v is the measure in (1.16) 
Moreover bi{n) = n~‘^, 62 (’^) = and = c; 


D is given by (1.15) and c,? = 1 


b) otherwise, under the assumptions (1.9) and (l.lf), hi{n) = (r(n)) , 62 ( 0 ) = {r{n)p{n)) 


Proof. To simplify the notation we just write hi and 62 - We prove the convergence of Laplace 
transform ( 61 T 0 , b 2 Mi) under the measure First, recalling the dehnition of the conditional 


measure given in (3.1), an elementary calculation using the branching property implies 


(4.27) 
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Thanks to Lemma 4 of [3] and Lemma 4.6 


(exp {-sbiTo - tb 2 Mi)) -)■ exp (-x i (1 - e ) , 


a(n) 


(4.28) 


so it remains to calculate the limit of the second factor in (4.27). Using again [3] together with 
the equality (4.17), we have 


(exp (—sfeiTo — tb 2 Mi)) = 1 — st 


e ^yftnib^^x,b 2 ^y)dxdy, 


'0 JO 


as before iJn{x,y) := Pi*'’^^(To > x or Mi > y). This implies due to Lemma 4 (ii) in and 
Lemma 4.6 (ii) 


(exp (—s6iTo — tb 2 Mi) Mi) - )-s 


'0 JO 


e ( X A — ) dxdy 


— st 


'0 JO 


ye ^x A —^ dxdy. 


Computing the integrals we get 

(exp {—sbiTo — tb 2 Mi) Mi) 

This finishes the proof. 


^-{s+c^t)z 


z'd{dz) 


(4.29) 


We can now continue with the proof of Theorem 1.9, assuming that (4.26) holds for k and 


to verify the case k + 1. Let be Jhc = o'((Mj_i,Tj), j < k) and (Tq,M() an independent copy of 
(To, Ml). Recalling the dehnition of the measure we have. 


pt ^g-AoTo-0iMi ^ ^ ^ g-AfeTfc— 



£3 A/c —l^fc —1 _ \^P 



Then applying the identity (4.27[) we get 


pt AqXo—S iMi _ _ _ g—AfcTfc—j 

= PP ^e-^oTo-OiMi . . .g-Afe_iTfc_i-0feMfe^ j ^ ^g-AfeT'-efc+iMj j ^ _ 

Now using the Markov property of {(T„, Mn+i),n G Z+} and writing the terms suitably, we get 

pt ^g-AoXo-6»iMi _ _ _ g—Afcrj.-0j,i_iMfei_i j 

T^-dk+iM[ 

-(ek-b2logEP_^ Mk 


^g-MogEP(e 


X E^ I (i~^oXo-diMi _ _ _ g-Afc_iTs,_ig 
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Taking in the previous equality (Ai_i, 6i) = (6iSj_i, b 2 ti), i = I,k + 1 and (a,p) = (a(n),p(n)), 


we use hypothesis (1.13) and (1.14) to obtain 


-p(n)t 

~a{n) 




ETC I 


-biSi-iTi-i-b2tiMi 


^ Z =1 


^ [pPWl ^g-6iso7b-fe2tiMi _ _ _ g-bi^ffc-iTfe_i-fe2(tii+K(^ffc+ctfe+i))Arfe^ 

Now we have to calculate the limit of each factor. The hrst one converges towards to 1 thanks to 


(4.28). Besides, to get 


Cg-feis,T'-fe2t.+iM(^A -^ i{Ko{Sk + C^h+l)), 

\ / n^oo 


(4.30) 


we use the convergence (4.29) together the convention fi:o(s) = s. As in Proposition 4.5, in order 


to conclude we use the induction hypothesis with s'_i +^ < i < k defined in (4.23). 


A Proof of Lemma 4.1 


A consequence of (1.9) is that the measure dehned on [0, oo) by 


/i(a;) := / ZTi^{z)dz, x > 0 


(A.l) 


is such that x i—)■ /i(x) is RV^ ". Then from the Tauberian-Abelian Theorem (see [1]), its Laplace 


transform G RVq and 


p(x) 


r(3 — a) 


T (1/x), X —>■ oo. 


Observe that 


A^T(A) = / (1 — e — Ay)7r~^(dy) + / A 2 :(l — e ^^)7r~^(dz). 

Jo Jo 


Then 


A2£^(A) = e(i- - Ae+e"^^^) , A ^ 0. 

As consequence of the definition of the measure y and the approximations above, 

7f+ (1/A) ~ c„E (l - , A ^ 0, 

where Cq, = l/r(3 — a). Hence for all x > 0, 


E(l-e-^«+ - ATe-^«^) 


T(x/A) 


7r+ (x/A) 
Ca7f+(1/A) A^o' 


> c„x 


(A.2) 


(A.3) 


(A.4) 
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We set r{n) = 1 —e j and define the measnre on (0, oo), mn{dy) = 

r{n)Tt~^{ndy). The convergence in (A.4) implies 

, c„ dy 


mn{x,oo) 

Therefore, for all 0 < x < ?/ < cxd 


mn{x,y] 


a y 


-)■ 


l+a ’ 


for all X > 0. 


Cq, dz 

yl-\-a 


This implies that the measnre on (0, cx)) dehned by mn{dy) = r{n)n~^{ndy) converges vagnely 
towards We also have 

[ y‘^l{y<a;}r{n)Tr+{ndy) -^ c„ [ y^l{y<x} 


l+a 


y 

Using an argnment of monote class to dednce the above convergence over intervals I C (0, cxd). 
Thns we obtain the convergence of the Laplace transform of the measnre /r. This complete the 
proof becanse /i is regnlarly varying at inhnity with indice 2 —a and its Laplace transform satishes 


the identity (A.2). 


B Proof of Proposition 4.8 


Before proving Proposition |4.8[ we wonld like to present some basic aspect of fnnctional conver¬ 
gence of stochastic process, farther details can be fonnd in [S]. Is well know that the law of a 
Levy process {X^ : f > 0} on is determined by that of random variable Xi, which is inhnitely 
divisible random variable, and according to the Levy-Khintchine formula has characteristic expo¬ 
nent 

T(u) = iu ■ b - ■ cu^ + J - 1 - iu ■ h(x)) F(dx) 

where b G R'^, c is a d x d symmetric nonnegative matrix, F is a positive measure on R'^ with 
F({0}) = 0 and /(I A |xp)F(dx) < oo, h is a truncation function from R'^ to R'^, that is, bounded 
measurable satisfying 

h(x) = o(|x|), |x| —)■ 0. 

Hence an inhnitely divisible distribution, and therefore a Levy process, is uniquely characterized 
by the triple (b, c, F). Another useful related equality is a d x d symmetric nonnegative matrix, 
called the modihed second characteristic, and dehned as follows 

_|_ J /i*(x)h-^ (x)F(dx), i,j = 1,2, ...,d. 

According to Theorem VIL2.9 of [8], if {Fn}n>i is a sequence of inhnitely divisible distributions 
on R'^. Then F^ ^ F weakly if and only if 

bn b 

Cn ^ C 

Fn(g) F(g) for all g G C'i(R'^), 
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where C'i([R'^) is a convergence-determing class for the weak convergence indnced by all continnons 
bonnded non-negative fnnctions —)■ R, vanishing at the origin and with limit at infinity. We 
stress that here we consider characteristics relatives to a continuous truncation function, h. 

In a more general sense, a d-dimensional semimartingale W, has associated a characteristic 
triplet (B, C, z/) consisting in: 

- B = {B^)i<d a predictable process with components of finite variation over each interval 

| 0 ,(|. 


- C = a continuous process, namely 

(Jij = 


where is the continuous martingale part of W. 
- V a. predictable random measure on R+ x R*^. 

A second modified characteristic C is also defined. 


+ (tfV) 

S<t 


h*(x)z/({s} X dx) 


h^(w)h'{{s} X dw) 


If W has no fixed times of discontinuity, in which case B is continuous, and \h{x)\‘^ * vt < oo, it 
reduces to 

Ci^ = Ci^ + (h^h^) * ut. 

According to Theorem VII.3.4 of [S], the necessary and sufficient conditions to assure the func¬ 
tional convergence of a sequence of semimartingales Vh"' towards W are given also in terms of 
their characteristics: 


sup \Bg — Bs\ —t 0, for all f > 0, 

S<t 

& C, for all t e D, 

g * z/” -)■ g * z/i, for alH G A), g G C'i(R'^), (B.l) 


where D is a dense subset of R_ 


We now turn to prove the convergence claimed in Proposition 4.8 In order to apply the 


Theorem VII.3.4 of |8], hrst we will prove the convergence of the characteristics of the process 




(n)t 


(r(n)t/n,0), t>0, 


(B.2) 


St = ^ (d“V«.d”’‘V''(n)p(n)) , keN, 

i=l 

and {Nt,t > 0} is a Poisson process with parameter one, independent of the sequence 

The following Lemma establishes the previous statement. We will use this result as a device to 
study the characteristics of which are closely related to those of — (r(n)t/n,0). 



Lemma B.l. The process defined in (B.2) is a semimartingale with characteristics relatives to 


a continuous truncation function h given by 

b” = r(77,)tE [h — {r{n)t/n, 0), 

= 0, = r{n)tE [K h, J = 1, 2, 

Ffi[d^ = r(n)t7r((ix), 


(B.3) 


where b{n)f^"''> = /n, r{n)p{n)'j, 7r*^”'^((ix) = G G dx 2 )- Moreover, 

in the regime and ( fi. we have the following weak convergence in the sense of Skorohod 

topology 

(B.4) 


'N, 


r(n)t 


(r(n)t/n,0) :t>0] ^ (Xt,t : t > 0), 


where Xt is a spectrally positive a-stable process with parameter a G (1,2). In particular, we 


obtain the convergence of the characteristics in (B.3) towards those relatives to {{Xt,t) : t > 0) 
and characteristic exponent Ca|A|", that is 


hf = [t 


'(O,oo) 


HHv) - y)cay ), 


Ff = 0, c*/ = E [h, (Xi) h, (Xi)] j = 1,2, 

Ffidx) = tCaXi dxi5Q{dx2), 
where h is a continuous truncation function from P to P. 


(B.5) 


Proof of Lemma B.l Note that for u = (A,^) G P^, 


E 


iu-S" 




(B.6) 


Then the exponent in the righthand side of the previons eqnality can be written as follows 
til I u ■ h{b{n)x.)r{n)7i{dx.) + f f — 1 — iu ■ h(6(n)x)) r(n)7r*'"^((ix) 

(0,oo) J (0,oo) J (0,oo) J (0,oo) 

(B.7) 
is in- 


where 6(n)x = {xi/n,X 2 /r{n)p{n)),'K^'^\dx) = P(^h"') ^ G dx 2 ). From here 

hnitely divisible, also we can dednce that the characteristics of the process — {r{n)t/n, 0) : 


r(n)t 


t > 0} are given by (B.3). Thanks to Theorem II.3.11 of [H] this process is a Levy process and a 
semimartingale. 


Besides, to get the convergence in (B.4) we shall prove the convergence of the characteristic 


fnnctions. This fact is verihed using Corollary |4.3[ together with the fact that conditionally to 
^(+) _ ^ distribution of is Binomial with parameter (fc, 1—p(n)), as well as the assumption 
that has mean 1. 


'( 0 ,oo) V r{n)p{n) 

'y,, , y 


y. 


y 


p{n) J J r{n)Ti^{dy) 


+ i ( —(1 — p{n)) H— p(?7-) ] r{n)t 


Pi.n) ) - ( -(1 -p(n)) + ^ p{n) ) y ) r{n)7r^{dy) 

n r[n)p[n) 


n 


r{n)p{n) 
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Then by the continuity of h and the assumptions (1.13) and (1.10), the previous display becomes 


t [ /^e-((i-P(n))(i-eWn)_pH(i_ei«AWPW)), _ ^ fy\\r{n)7r^{dy) 

■/(0,oo) \n// 


+ i (-(1 - p(n)) + f M n)] r{n)t. 

\n r[n)p[n) J 

where h is the truncation function from IR to K obtained as projection of h in the second coordinate, 
that is, 

~ t j j'g-((l-p(n))(l-e‘Vn)_p(,)(i_ei«/Kn)p(u))),,j^ _ ^ r{n)TT+{ndy) 


' (0,oo) 


'(0,oo) 


+ ti { Xh{y) - { X{1 - p{n)) + 


9 


r{n)p{n) 


np{n) ) y ) r{n)T:^{ndy) 


+ i ( -(1 - p{n)) + ^ p{n) ) r{n)t. 

' n r[n)p[n) ' 


Finally we have the convergence 

^ / iU'(S" 

E I e 




-)■ e 




(B.8) 


(B.9) 


where Ca is a constant depending on a that appears in Lemma |4.1 Indeed, this is a consequence 
of the arguments in this Lemma together with the hypothesis (|1.13 ) 


qojoo) 


g-((i-p(n))(i-e‘Vn)_p(„)(i_peM~Mn))),,j^ _ ^ _ iXh{y)) r{n)Tr+{ndy) 


t / - 1 - iXh{y)) Cay 

^ J{0,oo) 


while the second adding converges to 


it / X{h{y) - y)cay 

J (0,oo) 

To end we observe that r{n) G RV^ with a G (1, 2] and p{n) ~ cn~^ implies that 

r{n)p{n) 


G 0. 


Then 


Tl n^oo 


i ( -(1 - pin)) + , ^ , . pM] r(n)f - i-r(?7,)f, 
\n r[n)p[n) J n 


converges to i6*f. From (B.9) the characteristic of {Xt,t) are given by (B.5). As a consequence of 
Theorem VIL2.9 of [8] we have the convergence of the characteristics. Finally the characteristic 
exponent of Xt is c^A". ■ 
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We have now all the elements to prove Proposition |4.8 


Proof of Proposition Thanks to Theorem 2.3.11 of |8], is a semimartingale with char¬ 


acteristics relatives to h : —)■ given by 


Bt = - ([r(n)tj/n, 0) 

^n,ij ^ ^ h, (6(n)e("))] - E [h, {b{n)^^-^)] E ) t,j = 1,2. 

g* jyf = [r(n)tj E (g for g > 0 Borel, 

No w on ly remains to verify the conditions (B.l) where the limit characteristics agree with these 
in (B.5). In this direction we recall that b” = r{n)tE [h — (r(n)f/n, 0) and observe 


b - b.l 


1-8^ — b.l < I [r(n)sj — r{n)s\ E [h -|- |(r(n)s/n, 0) — {lr{n)s\/n, 0)| -f |b 

Then nsing the properties of the floor function we obtain 

|B," - b,| < E [h (6(n)e^"))] + |(s/n,0)| + |b^ - b,|. 

Thus by the convergence of b” established in the previous lemma, together with r(n) —>■ oo 
sup|5,”-b,| <Erh + I(t/n,0)|-f |b"-bt| -^ 0, 

s<t 


hence we have the first condition in (B.l). In order to determine the second one, let 6”’* be the 
Tth coordinate of b”, i = 1,2. Once more, applying the properties of the floor function we have 


1 - 




r{n)t J r{n)t 


6r&r’'-^E [h2 


< + I ^ 


{r{n)ty 


n 


E [*2 (Kn)^'">)] . 


Also 


and 


1 - 


1 ^ \ ~n,22 ^ ^ /^n,22 ^ ~n,22 , ^ /■lti,2\2 

cp - ) <Cp <cp +TTT7W ^t ’ 

r[n)t j r[n)t [r[n)t)^ 

^ ^ pyn,ll ^ ~n,ll . ^ /7 n,l'.2 i i ^ 


-n,ll 

o - 


r{n)t J r{n)t 


{bry < , lyA bry+hr + -Mn) - 1 ]. 


{r{n)t) 




As consequence of the convergence b” —)■ bj, 6"’* converges for i = 1,2. Then the above inequalities 
impliy —)■ cf, because r{n) G RVff and a G (1, 2). It is easily prooved that also g*z/” —)■ g*Ft 

for all g using that Fp = r{n)t'n'{dx.) converges to Ft, as we proved in the previous lemma, together 
with properties of the floor function. ■ 
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